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^SJ ' Abstract 

Deformation quantization is a powerful tool to quantize some classical systems espe- 
i-S^ . cially in noncommutative space. In this work we first show that for a class of special 

Hamiltonian one can easily find relevant time evolution functions and Wigner functions, 
f~| . which are intrinsic important quantities in the deformation quantization theory. Then 

I based on this observation we investigate a two coupled harmonic oscillators system on 

the general noncommutative phase space by requiring both spatial and momentum co- 
ordinates do not commute each other. We derive all the Wigner functions and the 
corresponding energy spectra for this system, and consider several interesting special 
cases, which lead to some significant results. 
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O ■ 1 Introduction 



Deformation quantization [T], based on Wigner's quasi-distribution function [2] and Weyl's 
correspondence between quantum operators and ordinary c-number phase-space functions 
5^ \ |2]-[5], have be extensively studied (for recent reviews, see e.g. Ref. [6]-[9]). 

In recent years, there has been much interest in the study of physics in noncommutative 
space. The ideas of noncommutative space-time and field theories defined on such a struc- 
ture started already in 1947 |101 lllj. It came up again first in the 1980's, when Connes 
formulated the mathematically rigorous framework of noncommutative geometry |121 [13] . 
A noncommutative space-time first appeared in physics was in the string theory, namely in 
the quantization of open string [14] , and the noncommutativity of space-time also plays an 
important role in quantum gravity [15]. Also in condensed matter physics the concept of 
noncommutative space-time is applied, such as the integer quantum Hall effect [16]. Since 
the noncommutativity between spatial and time coordinates may lead to some problems 
with unitary and causality, usually only spatial noncommutativity is considered. So far 
quantum theory on the noncommutative space has been studied extensively, and the main 
approach is based on the Weyl-Moyal correspondence which amounts to replacing the usual 
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product by a *-product in the noncommutative space. Therefore, deformation quantization 
has special significance in the study of physical systems on the noncommutative space. 

Some works on harmonic oscillators in the noncommutative space from the point of view 
of deformation quantization have been reported [17]. Considering there are many physi- 
cal models based on coupled harmonic oscillators [18j-[23j. it is interesting to investigate 
the coupled harmonic oscillators on the noncommutative space. In Ref. [24] two coupled 
harmonic oscillators on noncommutative plane with space-space noncommutativity were 
studied, and some results were obtained. In the paper the authors changed the noncommu- 
tative problems into problems in the usual commutative space with a well-known coordinates 
transformation in phase space, and solved the problems with normal quantum mechanics 
method. But there will be some difficulties of diagonalizing the new Hamiltonian obtained 
by the coordinates transformation in phase space, and the results such as the Wigner func- 
tions are still beyond to the commutative phase space, had not been transformed to those 
of the noncommutative space. So in the present paper we will remain in noncommutative 
space and use a new method of deformation quantization to solve this problem. 

First we show that for a class of Hamiltonian with special form one can easily derive 
the relevant *-Exponential function which determines the time evolution of the systems, 
and all the Wigner functions may be obtained by a Fourier-Dirichlet expansion of the *- 
Exponential function. When using this observation to the coupled harmonic oscillators on 
a general noncommutative phase space with both the spatial and momentum coordinates 
being noncommutative, we derive the explicit form of the relevant *-Exponential function, 
and then get all the Wigner functions and the corresponding energy spectra of the two 
coupled harmonic oscillators system on the noncommutative phase space. 

This paper is organized as follow. In Sec. [2] we briefly review the fundamental concepts 
of deformation quantization of a classical system on the noncommutative space. In Sec. [3] 
we discuss a class of Hamiltonian with special form and present some practical formulas. In 
Sec. m we study the properties of a system consisting of two coupled harmonic oscillators on 
the noncommutative phase space, derive its energy spectra and all the Wigner functions. In 
Sec. [5] we discuss some special cases of harmonic oscillators system based on the results in 
Sec. 131 and compare the results with those we have known in the literature. The summary 
and concluding remarks are in Sec. [6] . 

2 Deformation quantization on noncommutative phase space 

Consider a 41) general noncommutative phase space, the coordinates of position and mo- 
mentum are denoted by x={xi,X2} and p={pi,P2} , and their corresponding quantum 
operators x and p satisfy the following commutation relations [25\ 

[xi,Xj] = leijH 1 [Vi^Pj] = ''^'^ij^ , [xi,Pj] = ^Sijfi , (1) 
where i, j = 1, 2 , and e is the antisymmetric matrix 

In the deformation quantization theory of a classical system in the noncommutative 
space, one treats {x,p) and their functions as classical quantities, but replaces the ordinary 
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product between these functions by the following generalized *-product [8] 



*h *fi *u 



exp <j — ( dp^ - dp^ dx, ) + —eij d^, d^^ + —eij dp^ dp^ \ , (2) 



here we have used the Einstein summation convention and also in the latter part of this 
article without additional indication. The variables Xj, pi on the noncommutative phase 
space satisfy the following commutation relations similar to ([T]) 

[^jj^^j]^, = i^ijH , [PiiPj]* ~ '^eijV , — i6ijh , (3) 

where the Moyal bracket is defined as [f, g]^ = f * g — g * f . 

The time evolution function for a time-independent Hamiltonian H of the system is 
denoted by the *-Exponential function, which is the solution of the following equation 

^Tt--:^ i/i ifi-;^ iz^ \-r^ f Ht\ 

" ^ Y ^ -^^ikOxk ^ Pj-^ + yejWPi jExp ( ^ 1 • (4) 

where j, k,l = 1, 2. Eq. ^ corresponds to the time-dependent Schrodinger equation, and 
the *-Exponential function can be expressed by [8] 



*r, (5) 



n=0 

where 

(i?*)" = H*H*---*H . 

n times 

The generalized *~genvalue equation is 

H*Wn = Wn*H = EnWn , (6) 

where W is the Wigner function and E is the corresponding energy eigenvalue of the system. 
Eq. ([6]) corresponds to the time-independent Schrodinger equation. The relation between 
the Wigner functions and the *-Exponential function is: 

Exp(^) = f; e-^^-^/'^W. , (7) 

^ ^ ^ n=0 

this expression is also called Fourier-Dirichlet expansion for the time-evolution function. 



3 A class of Hamiltonian with special form 

Generally in order to obtain the *-Exponential function or the Wigner functions for some 
systems, one may try to solve the differential equation ([1]) or ([6]) which contains five or four 
variables. In many cases it is very difficult, but when the Hamiltonian is some special form, 
the problem can be simplified. 
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Now let us consider a Hamiltonian which can be written as a sum of two perfect square 
parts 

(8) 



H = {a ■ X + h ■ pf + {c - X + d- pf , 



where a={ai,a2}, b={bi,b2}, c={ci,C2} and d={di,d2}, and the coefficients Oj, bi, Ci and 
di are arbitrary real constants. The time evolution equation @ with the Hamiltonian ^ 
will become 



d 



dt 



Ht 



in— Exp l — ]=H* Exp — 



Ht 



+ \ cAx. 



2 ' 2 



i^i- 



(9) 



After some straightforward algebras we arrive at the following form 



in — Exp 
dt 



^\ = (^H-k%-eHd^)Exp 



Ht 



(10) 



where H is the Hamiltonian ([8|), and 

k = (aidi + a2d2 - 6iCi - 62C2) ^ + (aiC2 - a2Ci) /i + (61^2 - 62^1) 
= {a ■ d — b ■ c) h + {a A c) jj. + {b A d) u . 

Analogously, the generalized *-genvalue equation ([6]) becomes 

H *Wn=V\^n*H = (h - k^d^ - k^Hd'^\Wn = EnWn ■ 



(11) 



(12) 



Obviously, Eqs. (jlOp and (jl2p will work for any value of the parameters /x and v, especially 
for /i = and 1/ = 0, which correspond to the ordinary quantum theory. 

Now the Eqs. (jlOp and (jl2p are much simpler than the original form ^ and ([6]), they 
only contain one variable H in the right hand side of the equations. The solution of Eq. 
(UnD is [8J : 



Exp 



Ht\ 



1 



■ exp 



H kt 



ih J COS (kt/h) 
Using the generating function for the Laguerre polynomials 



^ ' 71=0 



(13) 



(14) 



where Ln{z) is the Laguerre polynomial, with s = e ^i^*/^ and z = 2H/k , one can write 
the *-Exponential function (|13p as 



Expf:^") =f;2(-l)"e-'=*(2n+i)/n^-HA^^j'2^\ _ 

^ ^ ^ n=0 ^ ^ 



(15) 
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Comparing this expression with the Fourier-Dirichlet expansion for the time evolution func- 
tion ([7j), we get the Wigner functions and the corresponding energy eigenvalue^: 

W„ = e-^/'=L„ (^^^ , En = {2n + l)k. (16) 

They are also the solutions of Eq. ()12p . 

Obviously, for most physical systems their Hamiltonian functions are not exactly the 
form ([8]), but we can use it to simplify our calculation in some cases, for example in the 
case of coupled harmonic oscillators. 



4 Coupled harmonic oscillators on noncommutative space 

Now consider the two coupled harmonic oscillators system [21^ [24j on the phase space, the 
Hamiltonian can be written as 

^0 = A' + i^Pl + \ {CiXl + C2XI + C3X1X2) , (17) 

where mi, are masses, and Ci,C2,C3 are constant parameters. After rescaling the 
coordinates of the phase space 



-2;^ \m,) ^ ^^^^ 



one can rewrite Hn as 



-f^i = ^ [pi +P2) + \ {(^i^l + (^2x1 + 0^x1x2) , (19) 



where m, ci, C2 and C3 are 



m = y/mim2, ci = Cu/ — , C2 = C2W— , C3 = C3 . (20) 

V mi y 1712 



To remove the interaction term in Hi (jl9p . one may use the following transformation 

yi\ /cosf -sinfWxA AiA /^cos f -sinfWpA 
2/2; Vsinf cosf^VW' W cosf^VW' 

which is a unitary rotation with the mixing angle a on the phase space. Applying the 
transformation pip to ()19p . and letting the parameter a satisfy the condition 

tana = — , (22) 
C2 - ci 

one can find the Hamiltonian (1191) will be 



H2 = ^i<l! + gl) + ^{e%! + e-^V.), (23) 



^In this work we will always ignore the normalized constant of any Wigner functions. 
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where 

1 /Zr^ J „2r, _ Ci + C2 + V(C1 - C2)2 + 



^= oV^^i^2-c^ , e'' = , (24) 

and the condition 4ciC2 > c| must be fulfilled. Since the transformation ()2ip does not 
change the structure of the generalized *-product ([2]), the coordinates yi, 7/2, 91 and ^2 
should satisfy the commutation relations ([3]) the same as those of X\ , X2 , Pi and P2 ■ 
The Hamiltonian H2 ()23p can be separated into the following two parts 

/ e^\/^sina cos a \ / e~'^\/l{ sinb cos 5 \ , , 

/ e''\/K'cosa sin a \ / e"'^'^/K cosh sin 6 \ 
^2 = 7= yi 7= 92 + 7= 2/2 7= 91 • (26) 



^/2 ''^ V^^'j ' \ V2 ^ ^ 

Obviously, with any value of a and 6, there should be Jix + 7^2 = -f^2- 

Using the observation obtained in Sec. [3] , we will immediately get the expressions of 
the time-evolution functions for Hx ([25]) and "^2 

^"p^ (it J = ^^^(MM U ^ J ' ^''^ 

^"P^ (^) = cos(fc2t/?.) ¥) ' 

where 



(29) 



The Wigner functions and the corresponding energy spectra are 

= e-'^'/'^Lr,, , < = (2ni + l)h ; (30) 

= e-^^/'^^L^, (^) , < = (2n2 + 1)^2 • (31) 
It is easy to verify when a and b take the following values 



ki = — = I e ' smacos — e ' sm ocos a) H smasm b cos a cos , 

2x/m 2 2m 



Ko = — 1= e ' sm a cos — e' sm cos a H cos a cos sm a sm 

2^M ^ ^ 2 2m 



sin(a — 0) = , cos(a — 0) — 



. , Km(e^ — e ^) / ,n Kmii + v 
sin(a + 0) = , cos(a + 0) = 

P2 P2 



or 



1 / n^/K^ie^ - e"'') fiy/K^ie'^ + e"^) 

a = — arctan -— h arctan ■ 



2 I —{Km^ + u) Kmn — v 

arctan ; ^ arctan ■ 



2 I —{Kmfi + h') Kvajx — V 
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(32) 



(33) 



then Til (|25p will be commutative with TI2 (|26p under the Moyal bracket 

[Wl,W2]*=Wl* 7^2-^2*^1 = 0, (34) 
and their generalized *-product is equal to their ordinary product 

ni*n2 = W1W2 = n2*ni . (35) 

f3i and (32 in ([32]) are 



/32 = (e'i - e-^iyn^Km + (/Cm// + zy)2 , 



(36) 



or 



with 



= {e^i + e-'')?i/^Vl + Ai , /32 = (e'' - e-'')n/^v'l + ^2 , (37) 



here Ai and A2 denote the effect of the noncommutativity of the phase space. When /i = 
and = , Ai = A2 = , it returns to the ordinary commutative phase space. 
With (j29p and (j32p . ki and ^2 can be written as 



(39) 



ki = -^{f3i + /52) = ^ ((e" + e-'')^/r+Ar + (e^ - e-'')x/r+A^ 
fe2 = - /52) = ^ ((e" + e-'')7r+Ar - (e^ - e"") ^/T+A^ 



where oj = \J K/m . Then from the definition of the *-Exponential function ([5]) and Eq. (j35p , 
we obtain the time evolution function for the coupled harmonic oscillators H2 (I23p on the 
noncommutative phase space 

H2t\ f Tiit\ I Tl2t 



Expi ( ^ ) Exp '^^^ 



1 /Til ht\ 1 /W2, M\ 

exp I 7— tan -— — — — - exp I 7— tan -— , (4(Jj 



cos (kit/h) \iki h J cos{k2t/h) \ik2 h 
and the Wigner functions are 



the corresponding energies are 



En,n2 = E^n! + = i^^i + l)ki + (2n2 + 1)^2 

= 7r~('("'l + "-2 + 1)^1 + ("-1 - ?^2)/52 

2?Ti V 



2 ((ni + n2 + l)(e'' + e-'')yrTAT 



+ (m - n2)(e'' - e-") VT+A^ j . (42) 
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Through the inverse transform of (j2ip . all the results above can be easily expressed 
in terms of the original variables of the noncommutative phase space (xi, X2,pi,P2) or 
{Xi, X2, Pi, P2) ■ So, based on the observation in Sec. [3l we derive a satisfactory result 
of the two coupled harmonic oscillators system on the noncommutative phase space with 
deformation quantization method, obtain all the Wigner functions and the energy spectra. 



5 Some particular cases 

5.1 fi = and u = 

If = I' = 0, then Ai = A2 = 0, the noncommutative phase space will reduce to the 
ordinary phase space. In this case, one has 

(3^ = (e^ + e-'^)nVK^ , ^ = (e'' - e-'^)nVK^ , (43) 

and 



From ()32p we have 

sin(a — b) = 1 , cos(a — 6) = , sin(a + 6) = 1 , cos(a + 6) = , (45) 
so we may choose a = 7r/2 and 6 = 0. Then 7ii (|25|) and 7^2 (|26|) become 

and the Wigner functions are 

yu _ g-HiAi-W2/fe2^ C^\l 



exp |-^-^(e-''g? + e^ql) - ^ {e^yl + e-^yl) 



( 2 „ 2 '^VK^ _„ 2 \ /.^N 
X -_=e''g2 + — e "2/2 , (47) 



with the corresponding energy spactra 



En,n, = (2ni + l)h + (2n2 + 1)A:2 = (e''(2ni + 1) + e-'^(2n2 + 1) 

= nu:(e^{ni + \)+e~'^(n2 + \\] . (48) 



These results are exactly the same as those in Ref. 
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5.2 ^ = -v = e 

For simplicity, we also choose = K = m = 1 (e.g. let ci=C2=l and C3=0 in (j24p ). then 
and /ci, k2 become 



and 



sin(a — b) 



ki = k2 

n 



cos(a — b) 



(49) 



(50) 



but the value (a + 6) can not be fixed, so we may choose a = vr for simplicity, then we have 



1 fi hy2-eqiy^ 



+ q2 



1 (2 , (^91+^^2)^ 



(51) 



Thus the Wigner functions are 

W — p-'Hi/ki-n2/k2 T 



1 



exp 



X L 



2 



27^1 



yt + 



Lr. 



2^2 



(?i9i + M^ , {fiy2-0qi) 



ihy2 - Oq^f 



+ 



+ q2 



■III 



X L 



"2 



and the corresponding energies are 




(52) 



Er, 



^ ((ni + na + l)/3i + (ni - n2)/32) 
{ni + n2 + l)Vn^ + 9^ . 



(53) 



Therefore, in the case of ;U = — = 9, ()4ip and (I42p will reduce to the results of Ref. |17j . 
5.3 fi,u h 

When fj,, v <^ h , it leads to Ai, A2 <^ 1, and 



/ Ai , Ao 

Vl + Ai pa 1 + ^ , Vl + A2 « 1 + ^ 



(54) 



Thus the energy (j42p becomes 



nin2 



^ (ni + n2 + l)(e'' + e-'')v/l + Ai + (m - n2)(e'' - e"") VT+Aa 



(55) 



Comparing this result with the energy spectrum ()48p on the commutative phase space, we 
find that, if the noncommutativity parameters ^ and v are much smaller than the Planck 
constant, the noncommutativity of the space structure will cause a shift of the energy 
spectrum of the coupled harmonic oscillator system. So if we find a deviation of the spectra 
of the coupled harmonic oscillators system from the standard ones in some experiments, we 
may measure and determine the values of the noncommutativity parameters /i and v and 
investigate the effects of the noncommutativity. 

6 Conclusion 

In this paper, we first consider deformation quantization for a class of systems with special 
Hamiltonian. Using this method, we investigate the two coupled harmonic oscillators system 
on a general noncommutative phase space, and obtain the explicit expression of the time- 
evolution function for this system, and furthermore, derive all the Wigner functions and 
the corresponding energies. Since the study is for a very general noncommutative space, 
and the results are also rigorous and exact, so the results certainly can be reduced to 
specific expressions reported in the literature in different specific cases. In particular, the 
results show that when the spatial or the momentum coordinates do not commute with 
each other, the energy spectra of the coupled harmonic oscillators system will have a shift 
comparing to the standard ones, and this phenomenon indicates a possible way to explore 
the noncommutativity of the spatial or the momentum coordinates in future experiments. 
When both the spatial and the momentum coordinates are commutative with each other, 
the results return to those of the ordinary case. We also discuss a case in which the 
noncommutativity parameters \i and v are much smaller than the Planck constant, and 
calculate the quantity of the energy shift, so if we could compare the quantity with some 
relevant experimental data, we would get the bound of the noncommutativity parameters, 
which will be very significant. 

All these results come from our observation in Sec. [3l and we believe that this method 
should be useful in some other cases. 
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